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The scalar mode density perturbations in a the warm inflationary scenario are analysed with a view 
to predicting the amount of non-gaussianity produced by this scenario. The analysis assumes that 
the inflaton evolution is strongly damped by the radiation, with damping terms that are temperature 
independent. Entropy fluctuations during warm inflation play a crucial role in generating non- 
gaussianity and result in a distinctive signal which should be observable by the Planck satellite. 



t>. ! I. INTRODUCTION 

o 

. Observations of the cosmic microwave background are consistent with the existence of gaussian, weakly scale depen- 
dent, density perturbations as predicted by most inflationary models The amount of non-gaussianity produced by 
the simplest inflationary models is small and unlikely to be to be observable by the next generation of experiments, 
but this still leaves open the possibility that a slightly more exotic inflationary model could produce a measurable 
effect. 

One variation on inflation is the warm inflationary scenario (see also @). Warm inflation is characterised by 
the amount of radiation production during the inflationary era. If the radiation field is in a highly excited state 
during inflation, and this has a strong damping effect on the inflaton dynamics, then we have what is know as the 
strong regime of warm inflation. In this case, thermal fluctuations in the radiation are transfered to the inflaton 
0, IS HUH, 0] and become the primary source of density fluctuations. 

The possibility of warm inflation occurring in realistic particle models has been enhanced by the discovery of a 
decay mechanism which is present in many supersymmetric theories [lol [ill Il2| . In these models, the inflaton decays 
into light radiation fields through a heavy particle intermediary. If the coupling constants are sufficiently large, these 
£ — ' models can lead to warm inflation. 

Fluctuations in the cosmic microwave background allow us to measure the density fluctuations at the surface of last 
scattering. We know, in principle, how to evolve these fluctuations from early times using, for example, the Bardeen 
variable £ [13]. Observations can be compared to predictions for various moments of the probability distribution of 
£. The most important of these is the primordial power spectrum of fluctuations P((k), defined by the stochastic 
average 



o 



o 

§ ; (C(ki)C(k 2 )) - (2 7 r) 3 P c (fc 1 )^ 3 (ki + k 2 ). (1) 

^ . The bispectrum B^(ki, k 2 , fc 3 ), defined by 

(C(k 1 )C(k 2 )C(k 3 )) = (2 7 r) 3 i? c (fc 1 , k 2 , fc 3 )£ 3 (ki + k 2 + k 3 ), (2) 

can be used to examine the non-gaussianity in the density fluctuations The normalised amount of non-gaussianity in 
the bispectrum is described by a non-linearity function /nl, defined by 

, , , x 5 B(k 1 ,k 2 ,k 3 ) 

Jnl^MM) - e P( kl )P(k 2 ) + P(k 2 )P(k 3 ) + P(fe 3 )P(fci) ' [ ) 

where the 5/6 factor is convenient for cosmic microwave background comparisons (l4| . 

Non-gaussianity can arise during the inflationary era, and become a feature of the primordial density fluctuations. 
The amount of non-gaussianity produced in the simplest inflationary models is typically around a few per cent 
[HI . [lrj . [l7| . and can be related to a standard set of slow roll parameters [T^]. For comparison, the second order 
Sachs- Wolfe effect is expected to act as a source of non-gaussianity in the cosmic microwave background observations 
equivalent to f^L ~ 1 llSfl - 

Maldacena [201 ] introduced a simple argument which can be used to determine the non-linearity parameter for 
cold inflation in the 'squeezed triangle limit' k\ -C k 2 sa /c 3 . Density perturbations freeze out, i.e. their amplitude 
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becomes constant, when their wavelength exceeds the Hubble length. Perturbations with the smallest wave number 
k\ freeze out first. Their effect on the bispectrum is equivalent to rescaling of the wave numbers &2 an d k 3 in the 
power spectrum. The rescalling behaviour of the power spectrum is described by the spectral index of scalar density 
perturbations n, leading to the result 

f NL = ^{n-l). (4) 

This result is quite robust, and applies to many versions of inflation (2lj |. There are, however, reasons to be cautious 
when we try to a pply the same idea to warm inflationary models. In warm inflation, both the radiation and the 
inflaton fluctuate [22| . The non-linear coupling between these fluctuations acts as a source of non-gaussianity, and 
this source is impossible to describe in purely geometrical terms. Only one mode of fluctuation survives after horizon 
crossing, but by then the non-gaussianity is already imprinted on the curvature fluctuations. 

The non-gaussianity produced by warm inflationary models has been looked at previously by Gupta et al. (23l . l24j . 
Our approach builds upon such work, but now we include the nonlinear coupling between the radiation and inflaton 
fluctuations on sub-horizon scales and find that a far larger amount of non-gaussianity is produced. The contribution 
to the non-gaussianity found in the previous work was dependent on the third derivative of the inflaton potential, 
equivalent to a second order effect in the slow-roll parameter expansion. We shall see that there are large contributions 
to the non-gaussianity, appearing at zeroth order in the slow-roll approximation. 

The best observational limit on the non-linearity function at present is from the WMAP three-year data [lj, which 
gives l/Wil < 114. The Planck satellite observations have a predicted sensitivity limit of around |/atl| ~ 5 [lj. Our 
minimum prediction for |/at_l| lies well above the Planck threshold, with a distinctive angular dependence, and should 
provide a means to test warm inflation observationally. 

The paper is organised as follows. We begin in section II with a brief introduction to the notion of warm inflation. In 
section III, we introduce fluctuations of the inflaton field described by a Langevin equation. The density fluctuations 
on large scales, which affect observations on the cosmic microwave background are studied in section IV, followed by 
a discusion of some of the consequences of our results, and extensions, in the conclusion. Appendix A describes the 
first order perturbation theory relevant to warm inflation and appendix B evaluates some of the integrals encountered 
in the main text. 



II. WARM INFLATION 



Warm inflation occurs when there is a significant amount of particle production during the inflationary era. We 
shall assume that the particle interactions are strong enough to produce a thermal gas of radiation with temperature 
T. Warm inflation is said to occur when T is larger than the energy scale set by the expansion rate H. The production 
of radiation is associated with a damping effect on the inflaton, whose equation of motion becomes 

4> + (3H + T)0 + V4 = (5) 

where T((j),T) is a friction coefficient, H is the Hubble parameter and V$ is the <f> derivative of the inflaton potential 
V{4>, T). Energy transfer to the radiation is associated with an increase of radiation entropy density s r according to 
the equation 

Ts r + 3HTs r = T<j) 2 (6) 
The effectiveness of warm inflation can be parameterised by a parameter r, defined by 

When r> 1 the warm inflation is described as being in the strong regime. 

Temperature dependence in the friction coefficient is a feature of many models [Hj|, and can lead to interesting 
effects on the density perturbations [22j, but in order to simplify the account given here we shall make the following 
simplifications: 

• r = r(ct>) 
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In this case, the time evolution is described by the equations 

4> + ri> + v<t> = o, (8) 

p r + AHp r = Tj> 2 , (9) 

3H 2 = 4vrG (2V + 2p r + 2 ) (10) 



where p r is the radiation energy density. 

During inflation we apply a slow-roll approximation and drop the highest derivative terms in the equations of 
motion, 

r0 + y = o, (n) 

4Hp r = T<j) 2 , (12) 
3H 2 = 8nGV (13) 

The validity of the slow-roll approximation depends on the slow roll parameters defined in [22| , 

i fv^ 2 i fv^\ „ i fv^VA 



16ttG V V J ' V 8irG\V )' 13 8nG { TV ) (W) 

The slow-roll approximation holds when e <C r, r\ <C r and |?«r. Any quantity of order e/r will be described as 
being first order in the slow roll approximation. 

III. INFLATON FLUCTUATIONS 

Thermal fluctuations are the main source of density perturbations in warm inflation. Thermal noise is transfered to 
the inflaton field mostly on small scales. As the comoving wavelength of a perturbation expands, the thermal effects 
decrease until the fluctuation amplitude freezes out. This may occur when the wavelength of the fluctuation is still 
small in comparison with cosmological scales. 

The behaviour of a scalar field interacting with radiation can be analysed using the Schwinger-Keldysh approach 
to non-equilibrium field theory prl [27|. When the small-scale behaviour of the interactions is averaged out, a simple 
pict ure emerges in which the field can be described by a stochastic system evolving according to a Langevin equation 

- w 2 4>(x, t) + rfa, t) + = t), (15) 

where V 2 is the flat spacetime Laplacian and £ is a stochastic source. For a weakly interacting gas with T>r, the 
source term has a gaussian distribution with correlation function [29l | 

(£(x,t)Z(x',t')) =2TT5 {3 \x-x')6(t-t'). (16) 

We shall restrict ourselves to the gaussian noise source, although non-gaussianity in the noise could act as a source of 
non-gaussianity in the density fluctuations. 

We can use the equivalence principle to adapt the flat spacetime Langevin equation to an expanding universe with 
scale factor a. The rest frame of the fluid will have a non-zero 3— velocity v" with respect to the cosmological frame 
and we must include an advection term. The Langevin equation becomes 

4>(x, t) + T<j)(x, t) + r a - 2 v?d a <j>{x, t) + V^- a~ 2 d 2 4>{x, t) = £(z, t) (17) 

where d 2 is the Laplacian in an expanding frame with coordinates x a . The correlation function for the noise expressed 
in terms of the comoving cosmological coordinates becomes, 

(Z(x,t)£(x' ,t')} = 2TTa- 3 (2ir) 2 6 {3 \x - x')5(t - t'). (18) 

The inflaton will also generate metric inhomogeneities, but these can be small on small scales. In first order perturba- 
tion theory, with a suitable choice of gauge, the small scale metric fluctuations can be discarded on sub-horizon scales 
(see appendix We will therefore apply eq. (1 1 7[) on sub- horizon scales and use a matching argument to extend the 
fluctuations to large scales. 
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The stochastic equation for the inflaton field has some similarities to the one which is used in the theory of stochastic 
inflation [30| . and we shall adopt a method for analysing the fluctuations which resembles one originally used by 
Gangui et al. flU in that context. However, there are important differences, reflecting the different interpretation 
of the stochastic equation in the two applications. In stochastic inflation, the large-scale inflaton field is constructed 
only from modes which have wavelengths larger than the horizon, whereas the inflaton field used here is valid for 
wavelengths larger than the scale of interactions in the radiation, which is sub-horizon size. We therefore need to 
retain time and spatial derivative terms which are dropped in Stochastic inflation. 

We shall use a uniform expansion rate gauge. The analysis of the Langevin equation can be simplified by introducing 
a new time coordinate r = [aH)" 1 and using the slow roll approximation. We are led to the equation 

<f>"(x,r) - ^V(x,t) - 3ri£0 a 0(a;,T) + T^v^ - d 2 0(x, r) = (2TT) 1 / 2 £(x,t) (19) 

where a prime denotes a derivative with respect to r and we have kept only the leading terms in the slow roll 
approximation. The noise term has been rescaled so that its correlation function is now 

(i(x,r)i(x',r')) = {2-k)H^{x-x')8{t-t'). (20) 

This equation is non-linear because T, v r and T depend on (j>. 

Now we treat the source term as a small perturbation and expand the inflaton field 

(j>(x, t) = <P(t) + 5 1( j)(x, t) + 8 2 ct>(x, t) + ... (21) 

where 5\4> is the linear response due to the source £. This expansion is substituted into the langevin equation and 
then we take the Fourrier transform. Only the zeroth order terms in the slow roll approximation will be retained. 
The equations for the first two inflaton perturbations are 

L5x<p = {2TT) 1/2 £ (22) 
L5 2 4> = ^(2rT) 1 / 2 ( 5 1 T*e-3rr- 1 (fc tt Ur )*(fc Q 5 1 0) (23) 

where v r is the scalar velocity perturbation v" = —iak a v r , the operator L is defined by 

L 8 x <p = 5x(j)" - —5x4)' + k 2 8xi> (24) 

T 

and * denotes the convolution 

/*<?(k) = J " /(k-k^k'). (25) 

Terms involving SiT have not been included because they appear at first order in the slow roll expansion. Such 
terms produce non-gaussianity in the fluctuations which depends on the first order slow roll parameters. A useful 
consistency check can be performed to confirm that the terms which are first order in the slow roll expansion reproduce 
Maldacena's result (|4|) in the squeezed triangle limit. Gupta et al [23l l24j have analysed terms involving the third 
derivative of the potential V^^. These produce non-gaussianity which depends on second order slow roll parameters. 
The perturbation equations can be solved using green function techniques. The solution to Lu = j is 

1 r°° 

U=- G(kT,kT'){kT , ) 1 - 2U ]{kT')d T ' . (26) 

k J T 

The retarded green function G{z, z') can be found for v — T/(2H) constant in terms of Bessel functions, 

G{z,z')=-^z v z lv (J v (z)Y v (z')-J v {z')Y v {z)). for z < z'. (27) 



Corrections due to the time dependence of v are similar in size to terms which we have already discarded in the slow 
roll approximation. 
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A. Inflaton power spectrum 



The inflaton power spectrum P c j > (k,T) is defined by 

(ty(ki,r)ty(k2,r)) =P (fc 1 ,r)(2 7 r) 3 ( 5(k 1 +k 2 ) (28) 

Substituting the first order inflaton perturbation from (I22D . using the general solution (|26p and the correlation function 
HHJ) gives 

dT'G(kr, kr'fikr'f-^ 2T(t')T(t') (29) 

Integrals of this type are examined in appendix |5] There is a saddle point in the integral when v is large, which 
allows us to take the FT out of the integral and obtain 

P+(k, t) = fc~ 3 2T{t f ) T(t f ) F(k, t, t) (30) 

where r = Tp at the saddle point, or according to eq. (|B10j) . 

^3 



k = a S j"-{HTfl\ (31) 

We call the time in eq. (|31[) the freezeout time tp for the mode k. The freezeout time always precedes the horizon 
crossing time, which occurs when k = aH . 

The remaining integral F{k, r, r) is a special case of a more general expression 

/>OC 

F(k,T 1 ,r 2 ) = k dT l G{kT 1 ,kT , )G(kT 2 ,kT r )(kT r ) 2 - 4L ' (32) 

J To 



which is examined in appendix [Bj An analytic approximation of F valid for large values of v is given in eq. (|B5|) . 
With this we recover a result derived in [22| , 

r) = fc- 3 ^(iny/2 T f 1 + ^fcV + ..) (33) 



2 v ' V r 

Note that r decreases with time and P^(fc, r) approaches a constant value P<p(k) on a timescale set by the freezout 
time. Horizon crossing occurs at r = 1. 



B. Inflaton bispectrum 

The bispectrum of the inflaton fields is defined by 

(<ty(ki,T)<ty(k 2 ,T)<ty(k3,T)) = (2^) 3 B (fc 1 ,A: 2 ,fc3)'5 3 (ki+k 2 +k3). (34) 

The first order inflaton perturbations 8\<j) are gaussian fields and their bispectrum vanishes. The leading order 
contribution to the bispectrum must therefore include a contribution to 5<j) from the second order perturbation, 

(M(ki, r)M(k 2 , r)5 2 0(k 3 , r)) w (2 7 r) 3 J B (fc 1 , fc 2 , fc 3 )^(ki + k 2 + k 3 ), (35) 

cyclic 

where 'cyclic' denotes cyclic permutations of {ki, k 2 , k 3 }. The second order perturbation can be obtained by solving 
eq. (1231) . The most interesting effects arise from the thermal fluctuations which are responsible for the temperature 
perturbation 8\T and velocity perturbations v. 

There are two types of fluctuation in the radiation. The first type is the purely statistical type of fluctuation which 
is present due to the microscopic particle motions. These fluctuations have been analysed before [12]. They can 
be important, but they decrease rapidly after the freezout time defined in the previous section. The second type of 
fluctuation is driven by the energy flux from the inflaton field. These fluctuations are important for coupling the 
inflaton and radiation fields, and we shall consider these in more detail. 
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The fluctuations in the radiation field need only be evaluated to first order in perturbation theory. The first order 
perturbation equations for the complete system can be found in appendix [5J The momentum flux J from the decay 
of the inflaton field, given by 

J= -T4>5i(f>, (36) 

is particularly important. If we keep only the leading terms in the slow roll approximation, the energy density and 
velocity perturbations in a uniform expansion gauge satisfy the equations 

5xp r + mSxPr + (2OH 2 + ^k 2 a- 2 ^j 5i Pr = k 2 a~ 2 J (37) 

Sip r + 4H5ip r — — — fca -1 p r v r (38) 

o 

where the subscript k used in eq. (|A21[) to denote the gauge choice has been dropped. 
We can solve eq. ([37|) exactly using green function methods as before, 

S lPr = V3fcr 4 J™ r'~ 4 sin r) J{t') dr' (39) 

v r = - — fcr 4 (°° t'- a cos k(jl ~ T) J(t') dr 1 (40) 
4p r J T V3 

where t = (aH)^ 1 . We may obtain a reasonable approximation by taking the slowly varying terms outside the 
integrals. Using eq. (|3"6"j) and p r ocT 4 , we have 



8iT _ 3fcl> 
9I> f 

V r = / 

TPr Jt 

where, 



g r (kr,kT J )5i4>{T / )d7 J (41) 
g v {kT 1 kT l )5 1 <t>{T l )dT l (42) 



1 /z\ 4 z' — z 

9rM = - ? =(-) sin-^ (43) 



fhM = -s-^ (44) 

For later reference, we define the integral 

/°° z 5 f°° x 4 e~ zx 

g v (z,z')dz' = — J - + Sx 2 dx > ( 45 ) 

The fluctuations in the radiation dominate over all other terms which drive the second order inflaton fluctuations. 
We shall split the second order inflaton perturbation into a part 5 v <p driven by the fluid velocity and a part S r (j) driven 
by the thermal fluctuations, 

S 2 (j) = S v cj) + S r 4> (46) 

The associated parts of the bispectrum will be denoted by B^i and B r ±. We can substitute the thermal fluctuations 
(0U and |g2} into eq. (J23]) and use the slow roll eqs. (|1U13|) to get 



L 5 v (f> = 



A v t- 2 J dr' J ^ g v ( P T, P T>) p • (k - p) <^(p, t') M(k - p, r) (47) 



L<\<> = A, I «/r' / ^^SrOw,pT y )M(py)l(k-p,r) (48) 



where 



Ar = -~, A V =3(>1 (49) 

6 (0 (0 
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The three-point function is given by substituting the solution for S v <j) into eq. (|35p . 

(2 7 :) 3 B;(k 1 ,k 2 ,k 3 )S 3 (k 1 + k 2 + k 3 ) = E dT ' J X dT " J (Ss P ' (k 3 - P) x 

cyclic r r 

G(fc 3 T,A ; 3T')(fc3T') 1 ~ 2 M„(T')T'- 2 3t ,(pT', J5 T")fc 3 - 1 X 

(M(ki, r)(Si0(k 2 , r)«Si0(p, r")<^(k 3 - p, r')>. (50) 

The four-point function splits into a product of two-point functions. Using eq. (|45[) and the results from appendix 151 
we obtain an approximation valid for large 

Bl{k 1 MM)^Y, A,P^x)P^k2h^ r^'f^^ + ^Sl^^fcsrO^r') 1 ^. (51) 

cyclic 3 T \ \ J 2 } / 

From eq. (|B16[) . 

B%(k u k 2 ,k 3 ) » 18 ^L(r) V (P ( W(fc 2 )(fcr 2 + Afc 2 )ki • k 2 ) , (52) 

~ cyclic 

where r = 2v/3 = T/(3H) and an asymptotic series for the function L{r) is given by eq. (|B17[) . 
A similar analysis for B^ using eq. ()B5[) gives, 

BJ(fci,fe,fc 3 )»3| E ^(fcijP^fca) (-^2) 7 V (53) 

™ cyclic V \ 2 3 / / 

This approximation. is valid for large v. 



IV. DENSITY FLUCTUATIONS 



The small scale inflaton fluctuations freeze out well in advance of the time when they cross the horizon. Whilst 
these fluctuations are freezing out, the metric fluctuations are relatively small (in the uniform expansion rate gauge: 
see appendix |A"]) . This has two important consequences. In the first place, we are relieved of the arduous task 
of doing second order perturbation theory for the metric perturbations. Secondly, we do not have to consider the 
effects of metric perturbations in the thermal field theory used to obtain the Langevin equation (|17|) . By contrast, 
in Maldacena's analysis of non-gaussianity in cold inflation [2(3], it was necessary to consider the inflaton vacuum 
fluctuations in conjunction with metric perturbations because the second order perturbations where the same order 
in the slow roll parameters as the metric perturbations and arose on horizon scales. 

Eventually, the wavelength of the perturbations crosses the effective cosmological horizon the metric perturbations 
become important. On large scales it becomes possible to use a small-spatial-gradicnt expansion (first formalised by 
Salopeck and Bond [3l|). This approach allows us to define the curvature perturbation £ so that it is conserved even 
in the non- linear theory [32I l33l . l34j . The bispectrum and the non-linearity of the density fluctuations can be obtained, 
to a reasonable accuracy, by matching the small and large scale approximations at horizon crossing. 

The large scale behaviour is governed by the same equations as the homogeneous system. In particular, during 
inflation the slow roll equations can be used to relate the total pressure and density to the value of the inflaton field. 
In this situation, the fluctuations can be described entirely by the conserved expansion fluctuation £ on constant 
density hypersurfaces, which is defined for general hypersurfaces by 

C = iln(l + 2^)+ / (54) 
2 J p + p 

where tp is the spatial curvature perturbation. After using the slow roll equations, 

C = iln(l + 2v>) + J jd<t>, (55) 



where H = H(4>) and 4> = <fi(4>) are given from eqs. (| 1 lfl 13|) . 
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Consider a uniform curvature gauge cp = and £ = (((/))■ When the inflaton perturbations are expanded as before 
in eq. PTj) . we have 

C = Q>Sic/> + +C^2(j) + ^C^i^ *5i<j>+ ... (56) 

where 4> subscripts denote derivatives with respect to 4> an d C<f> = H/<p. Hence the power spectrum of density 
perturbations is 

P f (fe)=C^(fc) (57) 
where P<p{k) is evaluated in a uniform curvature gauge. The bispectrum, using eq. (|35|) . is given by 

B c (kx, k 2 , k 3 ) = C^(k u k 2 ,k 3 ) + ^(P (fc 1 )P (fc 2 ) + cyclic). (58) 

Note that, according to eqs. (|!Hll3p , C^^/C^ is first order in the slow roll expansion. This term is important for 
relating the non-gaussianity to the slow roll parameters, but in our case the non-gaussianity occurs at zeroth order in 
the slow roll expansion and we can neglect terms which are first order. 

In the previous section we evaluated the small scale inflaton perturbations in a uniform expansion-rate gauge. On 
sub-horizon scales, we have argued already that the metric perturbations are relatively small and we should expect 
that the inflaton perturbations in the uniform expansion-rate gauge and the uniform curvature gauge should be 
approximately equal. This can be checked explicitly at first order in perturbation theory using eq. (|A25[) in appendix 
[K\ We will therefore use our earlier results for the inflaton fluctuations when applying eq. (f58|) . 

The density fluctuation bispectrum can now be obtained by matching the large scale curvature bispectrum eq. (|58|) 
to the small scale inflaton bispectrum eqs. and (fB"5)) . The largest contribution to the bispectrum, which we 

denote by B^, comes from the velocity term B^, 

B v c (k u k 2 ,k 3 )^ 18 L(r) £ (P c (fc 1 )P c (fc 2 )(A; 1 - 2 + fc 2 - 2 )k 1 .k 2 ), (59) 

cyclic 

which is valid for large r = T/(3H), where L(r) is plotted in FigureQ] 

2-5 -i 



L 




r 



numerical 
asymptotic 
log fit 



FIG. 1: The function L(r) determines the size of the non-gaussianity. The upper curve is a numerical calculation and the lower 
an asymptotic approximation (|BlTf) valid for large r. A log-linear fit eq. (|B18f) is also shown. 



9 



The amount of non-gaussianity contained in this part of the bispectrum is described by the non-linearity function 
fiVL, defined as in eq. ([3]). If we take a scale free spectrum with P((k) oc k~ 3 , then 



t* (i, i ^ - 15 L ^ ( MiM + fc 2 2 )ki ■ k 2 fcf (fc 2 2 + fcg)k 2 ■ k 3 fc|(fc| + fc 2 )k 3 ■ ki 

JNl{K1, K 2> K 3) ~ p.p. p 1 pT2 + pp + pT2 

1 2 3 \ 12 2 3 3 1 



(60) 



The value of f^ L always lies in the range — 15L(r) < f^ L < (33/2)L(r). In the equilateral triangle limit, 

r NL (k,k,k) = -VoL(r). (61) 

The momentum dependence of the bispectrum contributes to the angular bispectrum in the cosmic microwave back- 
round. This would be an important signature of warm inflation. Figure [5] showns the momentum dependence of the 
non-linearity function f^ L on a contour plot. 




The derivation of our result may fail in the squeezed triangle limit, when k\ <C fc 2 ~ ks, because one of the modes 
may exit the horizon before the other modes have frozen out. The freezeout time is given by eq. (|3ip . so that the 
consistency condition on our result is that 



We should also be prepared for the possibility that the time evolution of r during the inflationary era gives an 
additional small dependence on k. 

The contribution from the radiation density fluctuation terms, which we denote by f^ L , is given by 



r NL (kuk 2 ,k 3 ) = -- - 3 ■ -3 ■ -3 12, + hrr-Si + ttvt2 (63) 




This is substantially smaller than f^ L , but has the advantage that it is independent of the parameters in the model. 
In the squeezed triangle limit, f^ L = —5/4, which is identical to the 'generic' result for the non-gaussianity produced 
by the curvaton model [34[ . 
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V. CONCLUSION 

We have presented a preliminary analysis of the amount of non-gaussianity in the density fluctuations in the warm 
inflationary scenario. The discussion has been restricted to the strong, thermal regime of warm inflation where the 
radiation produced from the inflaton vacuum energy thermalises and where there is a large friction coefficient T = T((/>) 
in the inflaton equation of motion. We have found that the interaction between inflaton and radiation fluctuations 
leads to a large non-gaussianity in the density fluctuations when compared to cold single field inflation. 

The amount of non-gaussianity in the bispectrum, measured by the non-linearity function Jnl with equal momenta, 
is given approximately by 

/ Ar£ «_15In(l + ^)-| J (64) 

where r = T/(3H) is the parameter which must be large for the strong type of warm inflation. There are no simple 
means of measuring the value of r from observations apart from this effect on the non-gaussianity. A naive use of 
the limit on the non-gaussianity parameter from the WMAP three- year data [l|, |/jvx| < 114, would correspond to 
r < 2.8 x 10 4 . 

The sensitivity of observations by the Planck satellite has been estimated to be around |/jvx| ~ 5, although this 
assumes that the primordial value of |/jvl| is independent of momentum [3]. A limit |/jvx| < 5, which might arise 
if Planck does not detect any primordial non-gaussianity, would correspond to r < 2.5 and would not be compatible 
with the strong version of warm inflation which we have assumed here. 

If warm inflation is realised in the way which we have supposed, then the prospects for observing primordial non- 
gaussianity in the cosmic microwave background are very good. The angular bispectrum of the cosmic microwave 
background could then be used to distinguish warm inflation from other causes of non-gaussianity, for example special 
cases of the curvaton scenario [34], [H| , which can also produce a significant amount of non-gaussianity. A useful way 
to parameterise the bispectrum would be to take 

lmax 

B c (k k , fa, fa) = ^2 J2 Pd k i) p d k 2)fi(ki, fci)P*(ki • k 2 ), (65) 

cylic 1= 

where Pi is a Legendre polynomial. Our result for warm inflation corresponds to the I = 1 term with fi linear 
in k. The curvaton scenario produces a non-linearity /o = $nl constant. If the angular dependence of the signal 
corresponded to the dominance of the I — 1 term, this would provide very strong evidence indeed in support of the 
warm inflation scenario. 

Further work can be done to generalise the results obtained here to more general types of warm inflation. One 
possible extension would be to consider temperature dependence in the friction term, i.e. T((f>,T). Numerical work 
in [22I ] shows that this increases the amount of interaction between the inflaton and temperature fluctuations, and 
therefore we might expect even more non-gaussianity to develop. So far, the behaviour of the perturbations in 
this situation is only understood numerically, and further analytic work would be desirable for predictions of the 
non-gaussianity. 

For a complete picture of the non-gaussianity in models of warm inflation, we should also consider the weak regime 
of warm inflation, where the friction coefficient is small. In the weak regime, the fluctuations freeze out at horizon 
crossing. This invalidates most of our analytic results, but it should still be possible to evaluate the necessary integrals 
numerically. 

APPENDIX A: FIRST ORDER COSMOLOGICAL PERTURBATIONS 

In this appendix we shall review the equations for the first order cosmological perturbations of an inflaton and 
radiation system. These where first constructed in refs. [36|, Our notation closely follows the review by Hwang 
and Noh [381 ]. 

The most general scalar perturbations of the metric can be written in the form 

ds 2 = -(1 + 2a)dt 2 - 2a 2 (3, a dtdx a + a 2 {5 a0 (l + 2<p) + 2^ aP )dx a dx p (Al) 

where a(t) is the scale factor and a, j3, tp and 7 depend on space and time. We always use a Fourier transform with 
wave vector k a to replace the dependence on x a . We will also find it useful to introduce the perturbed expansion k and 
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the shear x of the normals to the time slices. Their spatial Fourier transforms are related to the metric fluctuations 
by 

K = 3{-tp + Ha) + k 2 a~ 2 x (A2) 
X - a{f3 + aj) (A3) 

The energy momentum tensor splits into a radiation part T ra b, and an inflaton part T^, a b, given by 

Trab = (Pr + Pr)v ra V r b + Prdab (A4) 

T^ab = 4>,a<f>,b- {^9 cd (f>,c(f>.d + V)g ab (A5) 

where the p r , p r and v ra are the radiation density, pressure and 4— velocity. Energy transfer between the two 
components is described by a flux term Q a , 

Trab b = Qa (A6) 

The inflaton equation of motion given earlier is equivalent to the choice 

Q a = -Tv b r ct> ib (t>, a (A7) 

In the unperturbed system, the values of the total density p, pressure p and all other background quantities depend 
only on time. Perturbations to the density and velocity are defined by 

ST 6o = S P> ST 6& = ih&ip + P) v (A8) 

where the orthonormal basis is denoted by hats, and k a is the normalised wave vector. Perturbations to the radiation 
are defined in a similar way in terms of T rao . The perturbed density and pressure become 

dp = Sp r + <j)5(j) - (f) 2 a + V^Scj) (A9) 
Sp = wSp r + <j>5(j) - <j> 2 a - VfStj) (A10) 

where w — 1/3 and 6<p is the inflaton perturbation. Perturbations to the energy momentum transfer are described by 
the energy transfer SQ and momentum flux J, 

6Q 6 = -5Q, SQ & =ik & J (All) 

The explicit expressions obtained from eq. (|A7[) are 



5Q = 5T(j) 2 + 2T<j>5(j> - 2aY<j) 2 (A12) 

J = -T<j>5<j> (A13) 

The relevant first order Einstein equations become 

k 2 a- 2 Lp -Hk = 4ttG Sp (A14) 

K-k 2 a- 2 x = 12nGk- 1 a(p + p)v (A15) 

X + H X -a-Lp = (A16) 

k + 2Hk + (3H - k 2 a- 2 )a = AttG(Sp + 35p) (A17) 

The inflaton equation becomes 

54> + 3H5<j) + k 2 a- 2 5(j) + = <j)(n + a) + (20 + 3H<j))a - SQ (A18) 

The energy and momentum equations for the radiation become 

Sp r + 3H(1 + w)Sp r = ~ka~ 1 (l + w)v r + p r a + (l + w)p r n + SQ (A19) 

(a 4 p r v r )' = ka 3 p r a + ka 3 (l + wy 1 (wSp r - J) (A20) 

Only six of the seven equations (|A14IIA2"0|) are independent. There are seven functions (p, k, x, a, v r , 5p r and S(f>, but 
the gauge freedom t — > t + ((x) allows us to choose the functional form any one of them at will. 
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The choice of gauge is rather arbitrary, but for short-wavelength calculations where ak > H, the uniform expansion- 
rate gauge k — proves to be convenient. This gauge choice is indicated by a subscript k. If we eliminate the fluid 
velocity v rK , 

5p rK + (8 + iw)H5p rK + (3(1 + w)H + 15(1 + w)H 2 + wk 2 a^ 2 )Sp rK = 

k 2 a- 2 {J K - (1 + w)p r a K ) + a~ 5 (a 5 (<5Q K + p r a K ))' (A21) 

where a K is given by substituting eqs. (|A9|) and (|A10|) into eq. (|A17|) . 

{3H - k 2 a - 2 + 167rG(/) 2 )a K = 4ttG((1 + 3w)Sp rK + 4^4 - 2V <50 K ). (A22) 

In the uniform expansion-rate gauge, the metric perturbation a K drops out of eqs. (|A18|) and (|A21[) when we work to 
leading order in the slow roll approximation. This was implicitly assumed in the discussion of the stochastic inflaton 
equation in section III. 

Converting the perturbations to other gauges can be done by examining gauge invariant combinations. For example, 
the combination 

6<f> v = 6<f> - %<p (A23) 
ti 

is gauge invariant and represents the inflaton fluctuation in a constant curvature gauge. Hence, in terms of uniform 
expansion-rate quantities, 

H v = Scj> K - ^-tp K . (A24) 
ti 



We can substitute for tp K from eq. (|A14|) . and obtain the equation 

4irGa 2 <i 



Hk2 Sp K . (A25) 
Note that, during inflation, 5p K « V^SipK and we hnd that 8(j) v « Scj) K to leading order in the slow roll approximation. 

APPENDIX B: INTEGRALS 

We begin with an approximation to the integral 

/>oo 

F(k,n,T 2 ) = k dr'G(fcr 1 ,/cT , )G(fcr 2 ,fcr')(fcr') 2_4l/ (Bl) 

J T 2 

where v — T/2H and the retarded green function G is given in eq. (|27p . The leading terms for large v and fixed n, 
Ti come from 



F(k,n,T 2 ) « —k^T^kT^Y^kTl) / J v {kT') 2 T l2 - 2lJ dT' 

4 Jo 

This is a standard integral, 



(B2) 



n 



Ukr'Yr'^dT' = - v (B3) 

T R (v + l)T R (v - 1/2) 



where Tji is the gamma function. We also have 

7iTr(i/) V" ' 4^ 

Hence, 



*"n(*)~--^U + ^ + -")- (B4) 



7T 



ro ,„., J1 ..^i l+ ^ji> + ...,. 
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A slightly more chalenging problem is the integral 

/oo 
G(k lT , hT')G(k 2 T, hT'KhT'y-^ihT'y-^gir') &t> (B6) 

where g(f) is a smooth function. We proceed as above to get 

F(Ai,As 2 ,r)w^-(fciAs 2 ) 3/2 r 2, 'y w (fcir)y i ,(fc2T) / J v (hr') J lJ {k 2 r l )r a - 2v gir^dr 1 . (B7) 
4 Jo 

For large v there is a Debye approximation for the Bessel functions which is valid in the range < r < y{k\k 2 )~ x / 2 , 

J v {k lT ) ~ (27r^tanha 4 ) _1/2 e ,y(tanhQ ^ ai) (B8) 
where vsechcti = fcj-r' for i = 1, 2. The relevant part of the integral becomes 

F{k 1 ,k 2l T)K^{k l k 2 YT 2u Y u {k l T)Y v {k 2 T) V 2 - 2lJ / (cosha) 2 ^ 3 e^ tanhai - Ql+tanhQ2 - Q2 ) ff (r')da (B9) 
8 Jo 

where ^secha = {k\k 2 Y^ 2 T' '. We find that there is a saddle point in this range at the value of a corresponding to 
t = Tp j where 

T F = {<6u) 1 ' 2 (kl + k 2 2 )- 1 ' 2 (BIO) 

Expanding about the saddle point gives 

; ; \ 3 / 2 
k\k 2 • 



^-^)~V^Urt%J 9{tf) (B11) 

which agrees with our earlier result when g(r) = 1. This saddle point is responsible for the phenomenon of 'freezing 
out' of the thermal fluctuations. The value of r decreases with time and the fluctuations always freeze out before 
they cross the horizon at r = 1. 

One final integral which we require is 

/oo 
dT'G(hr, fc 1 r')(A;ir') 1 - 2 ^(fc2r'), (B12) 

where 

1 f°° x A e~ x 

9W = a/„ (B13) 

The leading order behaviour for large v is given by 

/"OO 

F v (k x ,k 2 ,r) ~ (k X T) v Y v (k x r)- \ k x dr' J v {k X T'){k X T') 1 - v g{k 2 T r ) (B14) 

* Jo 

We could use a saddle point approximation, and this produces a good approximation when g(z) is replaced by a 
constant. However, the integrand is too flat to produce a good enough approximation with the function g(z). Instead, 
we interchange the orders of integration and use the identity 



/>oo 

/ dz J l/ (z)z 1 - ,J (z 2 + 3x 2 ) 
Jo 



2\-l-ju _ 



iF 2 (l; 1 - /z, 1 + v;"&x 2 /A) 



2 1 +T fl (l + u)T R (l + v) 



+ 2^41(1 + , + lF2(1 + Ml 1 + M, 1 + - + M; 3x 2 /4) (B15) 
We can integrate over x, let [i — > and use the large ^ limit, to get the asymptotic expansion 

fci\ 2 1 /2^\ / fci \ 2 ( k\\ 1 



1 5) y - U £ v (B16) 
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where 

25 30 121 

L(r) ~ In 2r h7H ln2r h . . . (B17 

or r 

and 7 is Euler's constant 0.57721 . . .. The asymptotic expansion is not very accurate for moderate values of r, and a 
numerical evaluation is provided in fig [TJ We can use a logarithmic fit, 

L(r) « In (l+^) (B18) 

for rough estimates. 
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